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^ \ Abstract 

C ' Copositive linear Lyapunov functions are used along with dissipativity theory for stabiUty analysis and control 

of uncertain linear positive systems. Unlike usual results on linear systems, linear supply-rates are employed here 
for robustness and performance analysis using Li- and Loo-gains. Robust stability analysis is performed using 
Integral Linear Constraints (ILCs) for which several classes of uncertainties are discussed. The approach is then 
extended to robust stabilization and performance optimization. The obtained results are expressed in terms of robust 
linear programming problems that are equivalently turned into finite dimensional ones using Handelman's Theorem. 
Q ■ Several examples are provided for illustration. 
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O ■ I- Introduction 

(N 

Linear internally positive systems are a particular class of linear systems whose state takes only 
nonnegative values. Such models can represent many real world processes, from biology [1], passing 
^ through ecology and epidemiology [|2l, to networking flS). Compartmental models, used e.g. in biological, 
medical, epidemiological applications, are also generally expressed as (non)linear nonnegative systems 
IfH. Several works have been devoted to their analysis and control, see e.g. flU, [|5l, [|6l, [|71, [|8l, [|9l, [ITOl , 

im, ma, d. 

Quadratic Lyapunov functions of the form V{x) = x^Px, with P = P^ positive definite, are the 
most commonly used to study the stability of linear systems. Dissipativity theory with quadratic storage 
functions and quadratic supply-rates are also widely used for robustness analysis, e.g. through the full- 
block S-procedure [[T4l . and norms computation, e.g. the Hoo-norm and generalized 'H2-norm. This 
quadratic framework also allows to apply powerful analysis techniques such as those based on Integral 
Quadratic Constraints (IQCs) [[TSl (exploiting the KYP-Lemma llT6ll . ifTTll . fTSll and the Plancherel 
Theorem), and well-posedness theory [fT9ll , ll20l . [|2T1l which is mainly based on topological separation 
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|fT9l . An interesting point in these approaches is that many of the obtained resuks can be represented 
as optimization problems involving Linear Matrix Inequalities (LMIs) ll22l . a wide class of convex 
optimization problems solvable in polynomial-time using e.g. interior-point algorithms [|23l . 

Unlike general linear systems, the stability of linear positive systems can be losslessly analyzed by 
considering a positive definite diagonal matrix P. Such functions are then called diagonal Lyapunov 
functions and the corresponding stability notion diagonal stability [|24l . This particular form dramatically 
simplifies the structure of the stability and stabilization problems, allowing for instance the design of 
structured and decentralized controllers in a simple way (convex formulation). Convexity does not hold 
in the general linear systems case where the design of structured controllers is known to be a NP-hard 
problem [|25l . Recently, a KYP-Lemma for positive systems using diagonal Lyapunov functions has been 
obtained in ifTOl by using the results reported in [|9l. In the same vein, exact losslessness conditions 
for robust stability analysis have been further obtained in [fTTl where several classes of uncertainties are 
covered. 

There exists however another class of Lyapunov functions leading to necessary and sufficient conditions. 
They are referred to as linear copositive Lyapunov functions [|6l, [|71, [|26l , [|27l , [|28l , l|29ll and write 
V{x) = X^x, where A is a vector with positive entries. In such a case, the resulting stability condition 
can be expressed as a linear programming problem (convex again) and solved in an efficient way. Since 
the Lyapunov function is linear, there is no more relationship with the vector 2-norm and the L2-norm 
as in the quadratic case, but rather with the vector 1-norm and the Li-norm. This framework is then 
more suitable for the analysis of the Li-gain of positive systems and its consideration in robustness and 
performance characterization. 

In this paper, the stability analysis and control of uncertain positive systems is considered in the Li- 
induced norm and Loo-induced norm using linear copositive Lyapunov functions and dissipativity theory 
||30l , H with linear supply-rates. Stability analysis and control synthesis results for unperturbed systems 
are first provided to set up the ideas and introduce the important tools. It is shown that computing the Li- 
and Loo -gains for linear positive systems is tantamount to solving a linear programming problem, with 
a complexity growing linearly with respect to the system size. The Li-gain is determined via a direct 
application of dissipativity theory while the Loo-gain is computed as the Li-gain of the transposed of the 
original system. It is also shown that the computed gains are valid regardless of the sign of inputs and 
states, relaxing then the concept of positive system to certain conditions on the system matrices only. 
The consideration of nonnegative states and inputs are for computational and theoretical considerations 
only, i.e. the use of linear Lyapunov functions and linear supply-rates. Convex necessary and sufficient 
stabilization conditions using full, structured and bounded state-feedback controllers and accounting for 
performance bounds are then obtained. While the question of determining under which algebraic conditions 
a system can be made positive in closed-loop is still open, the proposed methodology can be applied to 
implicitly characterize them, together with a performance constraint. 

Robust stability analysis is performed using Linear Fractional Transformations (LFTs), a classical tool 
of robust analysis [1311 . very few used in the context of positive systems ll32l . [fTTIl . [fT2ll . The advantage 
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of using LFT is that parameter-varying systems depending rationally and polynomially in the parameters 
can be expressed as a simple interconnection of a time-invariant linear system and an uncertain matrix 
depending linearly on the parameters, simplifying then the analysis. The use of LFT also extends to any 
type of uncertainties and nonlinearities, e.g. delay operators, uncertain stable transfer functions, more 
general bounded uncertainties and static nonlinearities. In the context of uncertain positive systems, the 
overall system is rewritten as a positive interconnection of a nominal system and a matrix of uncertain 
positive operators. These uncertain operators are characterized through Integral Linear Constraints (ILCs), 
the linear counterpart of IQCs. Although the provided framework does not enjoy the availability of the 
KYP Lemma nor the Plancherel Theorem, a frequency domain analysis can still be used in order to 
select the scalings accurately. Several classes of uncertainties are discussed and it is shown that for linear 
time-invariant uncertainties, ILCs fully characterize their static-gain matrices. Based on this fact, several 
general robustness properties are discussed. 

Robust stability analysis results are finally derived using dissipativity theory and formulated as robust 
linear programming problems. Exact stability conditions are provided in the particular case of LTI positive 
uncertainties with fixed static-gain matrix. To improve tractability, Handelman's Theorem [|33l is used to 
produce linear programs involving a finite number of constraints. A procedure to reduce the number 
of extra variables introduced by Handelman's Theorem is also proposed to reduce the computational 
complexity of the approach. The results are finally extended to robust stabilization. It is shown that the 
presence of scalings does not destroy the convexity of the problem as it is often the case in full-block 
S-procedure [|T4l or IQC approaches. Several examples demonstrate the efficiency of the approach and its 
exactness for a certain class of uncertainties including delays. 

Outline: Section HI] introduces the problem, fundamental definitions and results. Section HlI] is devoted 
to the stability analysis of unperturbed systems and Section |IV] deals with their stabilization. This is 
extended to robust stability analysis and robust stabilization in Sections IVl and IVH respectively. In Section 
IVIIl Handelman's Theorem is used to relax the robust optimization problems formulated in the previous 
sections. The results are then finally illustrated through examples in Section IVIIIi 

Notations: 1„ G M" denotes the column vector containing entries equal to 1. For general real matrices 
or vectors A,B^ W^^^, the inequality A < {<)B is componentwise. Let x E M", the vector cr-norm, 
a positive integer, a < oo, is denoted by ||2;||cr = (X]r=i l^ir)^^*^' while the vector oo-norm is defined 
by Halloo = maxjg{i^..._„} \xi\. Given v : [0, oo) — M", the Lo--norm Ht'lli^, cr positive integer, a < oo, 
and the Loo-norm ||f||L^ are defined by \\v\\l^ = {J^ \ \v{t)\\'^dtY^'' and \ \v\\l^ = ess supi>o ||t^(t)||oo, 
respectively. The spaces of signals v : [0, oo) — > having finite Lo.-norm is denoted by L". When the 
dimension of the signal has no importance, we will use the shorthand L^. For a matrix X E jR^x™^ 
[X]r^i and [X]c^i denote the scalar entry, the i*'^ row and the i^^ column. A linear map i{x) = c^x 
is said to be copositive if (Fx > for all x > 0. We also define the sets = {a G M" : a > 0}, 
:= {a G : a> 0, ||a|| ^ 0} and := {a G : a > 0}. 
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II. Preliminaries 
Let us consider general LTI systems of the form: 

x{t) = Ax{t) + Bu{t) + Ew{t) 

z{t) = Cx{t) + Du{t) + Fw{t) (1) 

x(0) = Xq 

where x,xq E M", u E M™, w E W and z E are respectively the system state, the initial condition, 
the control input, the exogenous input and the controlled output. When the system ([T]) with w = is 
asymptotically stable, it defines an operator from L^- 3 w z E L^, a positive integer. Such an operator 
framework is suitable for defining, computing and optimizing norms of systems [1341 , and determining 
robustness and performance properties. 

Definition 1: Let us consider the uncontrolled version of system ([1), i.e. u = 0. The system is said to 
be positive if the following conditions hold: 

i) The matrix A is Metlzer, i.e. it has nonnegative off-diagonal entries, 

ii) The matrices E, C, F are nonnegative, i.e. they only have nonnegative entries. 

A consequence of the above definition is that a) when it; = 0, we have x{t) E M" for all t > and all 
xo E M" ; b) when xq = 0, we have z{t) E R\ for all t > and all w(t) E M+. These facts justify the 
denomination of positive system. Note that the terminology also degenerates to the case of autonomous 
positive systems for which only the matrix A must be Metzler. 

Definition 2: A linear map V{x) = X^x with V{Q) = is said to be a linear copositive Lyapunov 
function for the positive system x{t) = Ax{t) if both V{x) > and V{x) < hold for all x E M" . 

Definition 3 (L^-gains of operators): Given an operator S : i— )■ L^, a positive integer, the Lo.-gain 
||S|U^-L^ is defined as 

:= sup \\^w\\l,. 

Equivalently, it is the smallest 6 > such that 
holds for all w E L^j. 

Definition 4: The Li-gain of an asymptotically stable linear time-invariant system H with transfer 
function h{s) = C{sl — A)~'^E + F mapping p inputs to q outputs is given by [34|: 

\\H\\l,_l,= max <^ V / \hij{t)\dt\ (2) 

where hij{t) is the impulse response from input j to output i. In the same way, the Loo-gain is given by 




\\H\\l^-l^ = ^max^^ { / \h,,{t)\dt } . (3) 

The Li-gain quantifies the gain of the most infiuent input since the max is taken over the columns. In 
contrast, the Loo -gain of a system is the max taken over the rows and then characterizes the most sensitive 
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output. Note that in the SISO case, the 2 norms obviously coincide, and so do all the Lo--norms, as 
pointed out in [fTSll . Another important fact, needed later, is the correspondence between the Li-induced 
and Loo -induced norms using the notion of transposed system: 

Proposition 1: Let us consider a system H with transfer function h{s) = C{sl — A)~^E + F and 
its corresponding transposed system H* having transfer function h*{s) = E'^{sl — A^)C'^ + F'^ . The 
Loo-gain of H is related to the Li-gain of H* through the equality: 

1 1 1 1 Loo-Loo = (4) 

Proof: The proof follows from the definitions of the transposed system and the norms. ■ 

Proposition 2: Given an asymptotically stable positive system H with impulse response h{t) E W^^, 
t > 0, we have 

=max|lj^/i(0)} and | |if | U^.^^ = max |/J(0)lp} (5) 

where h is the Laplace transform of h. 

Proof: Since A is Metlzer and Hurwitz, then we have e^* > for all t > 0. It is then immediate to 
see that the impulse response h{t) = Ce^^E + F is nonnegative as well. From the asymptotic stability of 
the system and © we hence have 

s=0 

= m 

and we get the result for the Li-gain. The Loo -gain expression follows from Proposition [1] ■ 
The above results then show that the Li- and Loo -gains of positive systems are intimately related to each 
others, and to the static-gain matrix of the system. The formulas given in Proposition [2] also suggest that 
a theoretical analysis of the gains is possible. 

in. Stability and Performance Analysis of Unperturbed systems 

In this section, nonconservative stability and performance analysis criteria for unperturbed systems are 
derived. It is assumed throughout this section that the system ([T]) is positive and that the control input is 
identically zero, i.e. u = 0. Similar results can be found in [|35l . [|32l . [[T2l . 

A. Li-gain characterization and computation 

Lemma 1 (Li-gain characterization): Let us consider system ([T) and assume it is positive. Then, the 
following statements are equivalent: 

i) System ([T]) is asymptotically stable and the Li-gain of the transfer w ^-^ z is smaller than 7. 

ii) System ([T]) is asymptotically stable and 1^ /i(0) < 7IJ where h{s) is the transfer function associated 
to system ([T])- 

iii) System dU) is asymptotically stable and tq{F - CA'^E) < 7IJ. 
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iv) There exists A G 



such that the linear program 



7lJ 



l^F < 



(6a) 
(6b) 



is feasible. 

Proof: It is immediate that ii) and iii) are equivalent since /i(0) = F — CA~^E. We then prove that 
iv) and i) are equivalent and, finally, consider the equivalence of iii) and iv). 
Proof of iv)^i) 

The proof relies on dissipativity theory for nonnegative systems [[301, lISI. Let us consider the copositive 
linear storage function V{x) = X^x with A G M"_,_ and the supply rate s{w, z) = '^\\w\\i — \ \z\\i, 7 > 0. 
Then, according to dissipativity theory, if the functional 



^-[{x^w^z) = V{x{t)) — I s{w{rf), z{r]))dr]. 



(7) 

is decreasing along the trajectories solution of the system ©, then the system is dissipative with respect 
to the supply-rate s{w, z) and the Li-gain is bounded from above by 7. Moreover, since A > holds, 
then asymptotic stability of the system also follows. The derivative of "H along the trajectories of ([T) is 
given by 

U = X'x{t)-^llw{t) + l'lz{t) 



A^A + IJC 



X'E - 7I J + t^F 



x{t) 
w{t)_ 

if and only if the left factor 



Since the signals x{t) and w{t) are nonnegative, then "H is negative on 
is negative, or equivalently if the conditions © hold. 
Proof of i)^iv) 

The necessity comes from Theorems 5.1, 5.3 and 6.2 of [6J. This can also be viewed from the fact 
that the S-procedure is lossless for any number of constraints in the linear case [|36l . [[37l . An alternative 
proof given in [[T2l relies on Farkas' Lemma (which is actually an alternative way of seeing the linear 
S-procedure [[371 Remark 4]). 

Proof of iv)^iii) 

A similar proof is given in [12]. Assume 4) holds, then the matrix A is Metlzer and Hurwitz, implying 
in turn that A^^ is a nonpositive matrix, i.e. A"^ < 0. By right-multiplying inequality ([6a1 ) by A~^, we 
obtain A^ > -l^CA'^ and after substitution into ([60), we get the inequality - CA'^E) < 7!^ 

which coincides with the one of statement iii). 

Proof of iii)^iv) The proof relies on an explicit construction of a Lyapunov function such that 
asymptotic stability of the system implies the feasibility of the linear program of statement iv). Assume 
iii) holds, then there exist v G R'f , and e > such that A < and 



t^F + {eiy' -I'^CA 



.T 



-1- 



E < 7I 



T 



(8) 
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Noting that the temi inside parentheses is positive since u > and A ^ < 0, then we can let := 
ev'^ — \^CA~^. By right-multiplying this inequality by A, we obtain 

\^A + llC = ev^A (9) 

which is negative from assumptions v'^ A < and e > 0. Since the left-hand term of is identical to 
the one of (|6al ). the implication of the feasibility of (|6al ) is proved. Reorganizing now the terms of the 
inequality ([8]), we get X^E — 7!^ + l^F < 0, which is identical to (|6b1 ). The proof is complete. ■ 

Remark 1: It is interesting to point out that despite of being computed with the assumption of nonneg- 
ative input signals w E Li and nonnegative state values, the determined Li-gain is valid for any input 
signal in Li and any initial state xq G M". This is due to the fact that the Li-gain has altemative definition 
©J which depends on the nonnegativity of impulse response only. Nonnegative input and output signals 
are considered for theoretical and computational convenience only. 

Lemma [H can be used to compute the exact Li-gain of any asymptotically stable positive linear system 
as follows: 

Algorithm 1 (Computing the Li-gain): The gain coincides with the optimal value of the following 
linear programming problem: 

minA,^ 7 s.t. AG W^^, 7 > and ([SD hold. 

The complexity of the above linear programming problem is quite low since it is a linear program 
involving n + 1 decision variables and 2n + p + 1 constraints. The complexity then grows linearly with 
respect to the size of the system. 



B. Loo- gain characterization and computation 

The following result is the counterpart of Lemma \T\ 

Lemma 2 (Loo-gain characterization): Let us consider system ([T) that we assume to be positive. Then, 
the following statements are equivalent: 

i) System ([T]) is asymptotically stable and the Loo-gain of the transfer w ^-^ z is smaller than 7. 

ii) System © is asymptotically stable and h{0)lp < 7!^ where h{s) is the transfer function associated 
to system ([T]). 

iii) System © is asymptotically stable and (F - CA^^E)lp < 7!^. 

iv) There exists A G M"_,_ such that the linear program 

AX + Elp<0 (10a) 
CX - 7lg + Flp < (10b) 

is feasible. 

Proof: The proof relies on Proposition \T\ By substituting the matrices of the transposed system into 
the conditions © of Lemma [B we get the conditions (flOl ). The equivalence between the statements is 
immediate from Lemma \T\ ■ 
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Similarly as for the Li-gain, it is possible to compute the Lgo-gain through an optimization problem. 
Algorithm 2 (Computing the Loo-gain): The gain coincides with the optimal value of the following 
linear programming problem: 



Unlike the Li-gain case, the number of constraints is 2n + q + 1 while the number of variables remains 
the same. Again, the computational complexity grows linearly with respect to the size of the system. 

IV. Stabilization of unperturbed systems 

This section is devoted to the stabilization of unperturbed linear systems via state-feedback control laws 
of the form: 



where the controller gain K belongs to one of the following sets: 

1) The set of unconstrained controllers /C := M"^^". 

2) The set of structured controllers JCc := {K e M"''" : [K]ij = 0, E Sc}, where Sc is the set 
of indices corresponding to entries in the controller gain. 

3) The set of bounded controllers ICb := {K e M™''" : < K < K+}, where K' and are the 
lower and upper bounds on the controller gain K, respectively. 

Unlike general LTI systems, the design of structured and bounded controllers is not a NP-hard problem 
[|25l . [l38l . [f39l . By indeed using a diagonal Lyapunov function, necessary and sufficient conditions for the 
design structured controllers can be easily obtained. We show here that it is also the case in the current 
linear setting and that it readily extends to the stabilization of positive systems with guaranteed Loo-gain. 

Lemma 3 (Stabilization with K E IC): Let us consider the closed- loop system (fTI)- (fTTI ). where (H) is 
not necessarily a positive system, with transfer function hci{s, K) := (C + DK){sI — A — BK)~^E + F. 
Then, the following statements are equivalent: 

i) There exists a controller matrix K such that the closed-loop system (fTI)-(fTTI) is positive, asymptoti- 
cally stable and the Loo-gain of the transfer w ^ z h less than 7 > 0. 

ii) There exists a controller matrix K such that the closed-loop system (fTI)-(nTI) is positive, asymptoti- 
cally stable and verifies /ici(0,iC)lp < 7!^. 

iii) There exist A E and /ij E M™, i = 1, . . . , tj, such that the linear program 



minA,^ 7 s.t. A E ]R'[+, 7 > and ([lO]) hold. 



u{t) = Kxit) 



(11) 



n 




(12a) 



i=l 



n 




(12b) 



[^]tj \ + [B]r,ilJ-j > 0, 2, j = 1, . . . , n, j 
[C], . A, + [D]r^iHj > 0, 2 = 1, . . . , g, J = 1, . . . 



n 



(12d) 



(12c) 
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is feasible. In such a case, a suitable K is given by 



K 



(13) 



Proof: To prove the equivalence it is enough to show that statements 1) and 3) are equivalent. The 
rest follows from the proofs of Lemmas [H and |2] The closed-loop system is given by 



(14) 



(15) 



x{t) = {A + BK)x{t) + Ew{t), 
z{t) = {C + DK)x{t) + Fw{t). 

Substituting the closed-loop system into (flOl) yields 

{A + BK)X + Elp 1 ^ 

Noting that K\ = Yl^=i ^'il^]c,u it turns out that the change of variable /Xj = [/r]c,jAj linearizes the 
problem and yields (|12al )- (|12b| ). To ensure the positivity constraint on the closed-loop system we need to 
impose A + BK and C + DK to be Metzler and nonnegative, respectively. Using the same procedure as 
in Q, these constraints are captured by inequalities (|12c| )- (|12dl) . ■ 
The above theorem can hence be viewed as an extension of [7j where no performance criterion is 
considered. The following proposition discusses the cases K E JCc and K E JCb. 

Proposition 3 (Cases K E ICc and K E JCb): These cases can be easily handled by adding supplemen- 
tary constraints to the linear program of Lemma |3] 

• The design of a controller gain K E ICc ensured by considering the additional linear constraints 
fj-i ... /im = 0, for all (2, j) G 5c. 

J ij 

Prescribed bounds on the coefficients of the controller gain are imposed, i.e. K E JCbhy considering 
the additional linear constraints [K']r,iXi < fJ-i < [K+]r,iXi, i = 1, . . . , m, E M™^", < 

K+. 

Note that in both cases the necessity of the approach is preserved. The above results are thus noncon- 
servative. Other constraints like asymmetric bounds on the control input and the consideration of bounded 
states can also be easily considered [|71. 



V. Robust stability analysis and robust performance 

Let us focus now on uncertain linear systems subject to real parametric uncertainties 6 E S := [0, 1]^, 
iV > 0, of the form 

x{t) = As{S)x{t) + Bsi6)u{t) + Esi6)wi{t) 

zi{t) = Cs{6)x{t) + Ds{6)u{t) + Fs{6)w,{t) (16) 
x{0) = xo 

where x E M", xq E M", u E M™, wi E and zi E are the system state, the initial condition, the 
control input, the exogenous input and the performance output respectively. We assume in this section that 
the above system is positive, that is, for all 6 E 6, the matrix As{6) is Metzler, and Es{6),Cs{S), Fs{6) 
are nonnegative matrices. We also assume that the system matrices are continuous functions of 6. 
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While this section mainly focuses on positive systems subject to real parametric uncertainties, the 
proposed methodology basically applies to any type of positive interconnections involving, for instance, 
delays, positive infinite-dimensional operators, time-invariant and time-varying uncertain positive oper- 
ators, static sign-preserving nonlinearities, etc. A justification for focusing on parametric uncertainties 
lies in the particular property of positive systems that (as we shall see later in Section IV-BI) only the 
static-gain matrix is critical for evaluating stability of interconnections. This property implies that many 
problems involving uncertain positive transfer functions equivalently reduce to problems involving constant 
parametric uncertainties. 

A. Positive Linear Fractional Transformation 
Using LFT, the system (fT6l) is rewritten as 



where the loops signals ZQ{t),WQ{t) G M"" have been added. It is important to keep in mind that we 
are considering robustness in the Li-norm, the tractability of which relying on the nonnegativity of the 
loop signals wq and zq. Representation ([17] ) is not unique, even when minimal, so the arising question 
concems the constructability of a so-called positive Linear Fractional Representation(hFR), i.e. an LFR of 
the form (flTT ) with nonnegative loop signals and positive operators, corresponding to the positive system 
(fT6l) . Positive interconnections are recurrent in the analysis of positive systems and have already been 
considered in several works 0], [1401 . [[T2l . [|32l . By picking suitable matrices for the LFR, it is always 
possible to make the loop signals nonnegative. 

B. Handling positive uncertainties via Integral Linear Constraints 

In robust stability analysis theories addressing two interconnected systems, stability conditions very 
often consist of two separate conditions: one for each subsystem. The first one, generally very precise, 
is used to characterize the nominal system which usually enjoys nice properties like linearity, time- 
invariance, etc. The second condition, often more difficult to derive when high precision is sought, is 
used to characterize the uncertain part. Several approaches have been developed to study the stability of 
interconnections like small-gain results [|3T| and generalizations [[T4| . well-posedness/quadratic separation 
11411 . En (based on topological separation arguments [|T9l ) and IQCs [[T5|| . IQCs are very powerful objects 
that are able to implicitly describe, with high accuracy, uncertain operators through the characterization 
of their input/output signals. They unfortunately do not fit in the current framework since they are based 
on quadratic forms, while the proposed approach considers linear ones. Inspired from this idea. Integral 
Linear Constraints (ILCs) are considered in this paper. 



Wo{t) 



Ax{t) + EoWo{t) + Eiwiit) 
Cox{t) + Foowoit) + Foiwiit) 
Cixit) + Fiowoit) + Fiiwiit) 



(17) 
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1 ) Definitions of ILCs and preliminary results: 

Let us consider first an uncertain operator S : Li i-t- Li belonging to a known set S. The main idea is to 
implicitly characterize the set E using ILCs ai] 

S g |s : (I^WI + ^lMt)\) dt>{), z = T,w, 1 = 1,..^. (18) 

The vectors ipi^i and Lp2,i are called scalings and must be chosen according to the set S. Note that the 
absolute values disappear when considering a set of positive operators and nonnegative input signals. 

One important feature of IQCs lies in the fact that, by virtue of the Plancherel Theorem, it is possible to 
express the inequality in the frequency domain, domain in which the tuning of the scalings may be easier 
and/or more precise. Finally, by virtue of the Kalman-Yakubovich-Popov Lemma and the S-procedure, the 
frequency domain conditions are converted back to time-domain conditions, taking the form of tractable 
LMLproblems. The main concern is that the Plancherel Theorem does not exist in Li and we cannot 
expect to switch from time to frequency domain in the same spirit as in L2. It is however still possible 
to consider the frequency domain as stated in the following result: 

Lemma 4: The following equivalent statements hold: 

i) The ILC 

r+00 

/ {vlw{t) + '^lz{t)) dt>0 (19) 
Jo 

holds for every pairs of positive signals (w, z) such that z = Sw, where S is a positive operator. 

ii) The algebraic inequality 

^'^w{0) + ip^ziO) > (20) 

holds where w and 2 are the Laplace transform of the signals w and z, respectively. 
Proof: The proof follows from noting that 

j;-{vJwit) + ^lzit))dt = J^'--[{^Jwit) + ^lz{t))e-^']dt 

= ipJwiO) + iplziO)- 



s=o (21) 



Assuming that the pair of nonnegative signals {w, z) is related by the equality z = 'E'^rpjW where S^^j 
is a linear time-invariant positive operator, we have the following corollary: 
Corollary 1: The following equivalent statements hold: 

i) The ILC 

P+OO 

/ {vlw{t) + ip2z{t)) dt>0 (22) 
Jo 

holds for every pairs of nonnegative signals (w, z), z = Tj\j,^w, where S^^^^ is a linear time-invariant 
positive operator. 

ii) The algebraic inequality 

i^l + (/^^S+„(0) > (23) 

'the absolute value is componentwise. 
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holds where T,j^j,j(s) is the transfer function corresponding to the linear time-invariant positive 
operator T.'^^j. 

Proof: Since S^^^ is a linear time-invariant positive operator, then the left-hand side of the inequality 
(|20l ) can be written as (pf + (/?^E(0) {^{O). The nonnegativity of the letter expression is equivalent to 
the nonnegativity of the term into brackets since 'u;(0) = w{t)dt > 0. This concludes the proof. ■ 
2) Generic robustness results and remarks: 

Lemma |4] and Corollary [1] show that, similarly as for IQCs, the scalings may be selected both in 
the time-domain and the frequency domain. In the LTI case, the scalings just have to be selected 
according to the static-gain matrix of the system, which shows that the problem reduces to a problem 
with constant parametric uncertainties since the static-gain matrix is a (possibly uncertain) real matrix. 
Another conclusion is that, since only the zero-frequency is important, there is thus no need for capturing 
the entire frequency domain by using frequency-dependent scalings as in the /i-analysis [|42|1. IQC-based 
techniques [[TSll or approaches using dynamic D-scalings [43 J. It is hence expected that a constant scaling 
matrix acting on the static-gain matrix should be enough for obtaining interesting results. 

Along these lines, it is proved in [ | 12 j that the use of the scaled Li-gain yields necessary and sufficient 
conditions for the characterization of any interconnection of LTI positive systems. Similar results are 
obtained in the L2 framework in [fTTI . These results can be understood through the fact that the stability 
of interconnections of LTI positive systems is equivalent to the stability of interconnections of their static- 
gain matrix, i.e. higher order dynamics have no impact on stability. The critical stability information is 
hence concentrated in the static part of the dynamical systems. 

Based on the above remarks, generic robustness properties are formalized below: 

Theorem 1: The following statements are equivalent: 

i) The linear time-invariant positive operator S satisfies the ILC (l22l) . 

ii) The static-gain matrix S(0) satisfies (|23l) . 

iii) The linear time-invariant positive operator E' satisfies the ILC (|22)) . where = + s0(s) 
and sQ{s) is any positive asymptotically stable proper transfer function. 

Proof: The equivalence between i) and ii) follows from Corollary [1] The proof of ii) =^ i) follows 
from choosing 6(s) = 0. To show the converse, it is enough to remark that S'(0) = S(0) and hence S' 
satisfies the ILC (l22l) for any positive asymptotically stable strictly proper transfer function 0(s). ■ 
The above result interestingly shows that the set of uncertainties that satisfies the ILC is typically very 
large since there is basically no limit on the magnitude of the coefficients of 0(s) acting on powers of s. 

C. Main results 

Based on the results and discussions of the previous sections, the main results on robust stability 
analysis of linear positive systems subject to uncertain parametric uncertainties can be stated. Both Li- and 
Loo-induced norms are considered. The conditions are expressed as robust linear optimization/feasibility 
problems that are difficult to solve directly. An exact solving scheme based on Handelman's Theorem is 
proposed in Section IVIII 



13 



Theorem 2: Assume there exist a vector A G M"^ and polynomials <^i{5), ^2{^) G MJ^° such that the 
robust linear program 

A^A + v9i(5fCo + lJCi < 
A^i?o + V^2(5r + ¥^i(5)^Foo + l^Fio < (24) 
A^i?i-7lJ + ¥^i(<^fFoi + l^i^ii < 

ifiiSf + ^2{5YA{5)>Q (25) 

is feasible for all 5 G <5. Then, the uncertain linear positive system (fT6] ) is asymptotically stable and the 
Li-gain of the transfer wi — )■ zi is smaller than 7 > 0. 

Proof: The proof follows exactly the same lines as for Lemma [T] Note however that the supply-rate 

s{w, z) = -^i{5fzo{t) - i^2{5fwo{t) + 7l>i(0 - l^^i(0 (26) 

has to be considered here. ■ 
It is important to mention that the above condition is unlikely a necessary condition since the linear 
Lyapunov function does not depend on the parameters. It is indeed well-known that quadratic stability 
(parameter independent Lyapunov function) is more conservative than robust stability (parameter dependent 
Lyapunov function). 

The following result addresses the case when the ILC can be saturated using constant scalings: 
Theorem 3: The following statements are equivalent: 

i) System (fTTI) with A = Aq, Aq > constant, is asymptotically stable and the Li-gain of the transfer 
w ^ z is smaller than 7. 

ii) System (fTTI) with A(s) LTI, asymptotically stable, positive and A(0) = Aq > constant, is 
asymptotically stable and the Li-gain of the transfer w ^ z is smaller than 7. 

iii) There exist vectors A G M"^, (pi,^p2 G M"" such that the linear program 

A^A + <^fCo + lJCi < 
A^Eo + + (^f Foo + l^Fio < (27) 
A^Fi - 7IJ + (^f Foi + l^Fn < 

+ ip^Ao = (28) 

is feasible. 

Proof: The equivalence between i) and ii) follows from Theorem [T] It is easy to show from dissipativ- 
ity theory that iii) implies i). This can also be proved by substituting (^J by — (/^I^Aq in the inequalities (l27l) 
and further substituting by — A^F'o(/ — AqFoo)^^ — 1^Fio(/ — AqFoo)^^. This leads to the inequalities 

A^(A + Fo(/-AoFoo)-iAoCo) + l^(Ci + Fio(/-AoFoo)-iAoCo) < 

A^(A + Fo(/-AoFoo)-iAoFoi)-7l^ + l^(Fn + Fio(/-AoFoo)-iAoFoi) < 

which are exactly the conditions for asymptotic stability and bounded Li-gain. By reversing the calcula- 
tions, we can show that i) implies iii). ■ 
This result can indeed be extended to the Loo-case, this is omitted for brevity. 
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Theorem 4: Assume there exist a vector A G M"^ and polynomials <^i{5), ^2{^) ^ IR*^" such that the 
robust linear program 

A\ + Cq^i{5) + Eitp < 
E];\ + ^2{5) + Fj^^,{5) + Fll.p < (30) 
CiA-7l, + Foi^i(5) + Fnlp < 

+ ^2{5f^{5) > 0, (5 G (5 (31) 

is feasible for all 5 G 5. Then, the uncertain linear positive system (fT6] ) is asymptotically stable and the 
Loo -gain of the transfer wi — t- zi is smaller than 7 > 0. 

Proof: The proof is based on the use the LFT system corresponding of the transposed system of (fT6] ) 
given by: 

x{t) = A^xit) + EoWo{t) + C]^wi{t) 
zo{t) = Cox{t) + Foowoit) + Foiwiit) 
z,{t) = Efx{t) + F,oWo{t) + F[,w,{t) 
Wo{t) = A{6)zo{t) 

where the matrices Fio and Eq are specific matrices of the transposed system. All the other matrices are 
those of systems (fT6l) and (fT71) . ■ 
Remark 2: It must be stressed here that the Linear Fractional Transformation does not commute with 
the operation of transposition. In other words, the transposed of an LFT system does not coincide, in 
general, with the LFT of the transposed system. Some matrices may indeed remain unchanged (non 
transposed) while some others are different. This has motivated the use of the 'bar' and 'tilde' notations 
in (|32|) . For instance, when the system depends polynomially on the parameters, we may have the equalities 
-^00 = -^00) Fqi = Fqi and Cq = Cq. 

VI. Robust Stabilization 

The robust stabilization problem with Lqo -performance is solved in this section. It is interesting to note 
that the presence of scalings never destroys the convexity of the approach as this may occur in some 
robust control approaches based on IQC's [[TSl or the full-block S-procedure [[T4l . 

As in Section |Wl it is not necessary that the open-loop be positive, the only requirement is the 
nonnegativity of disturbance input matrices, i.e. E(S),F(S) nonnegative for all S E d. The LFR of the 
closed- loop transposed system is given by: 

i(t) = A{K)S:{t) + So{K)wo{t) + £i{K)wi{t) 

Zo{t) = Cox{t) + J'ooWoit) + J'oiWii^t) ^^^^ 
Zi{t) = CiX(t) + J'ioU)o(t) + J'llWl(t) 

Wo{t) = MSfzoit) 
where A{K) = (A^ + B^Kf, £i{K) = (C? + D^Kf and 



£o{K) = {A' + B'Kf ... {A'^ + B'^Kf {Cl + D^Kf . . . (Cf + D'^Kf 



(34) 
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The integers ?7, 6* > are related to the system dependency on the parameters. 

Theorem 5: The controlled system (fT6])- (fTn ) with G /C is asymptotically stable if there exist vectors 
A G M"_f,, V5i(5), (/J2(5) G R"", /^j e M™, i = 1, . . . , n and a scalar 7 > such that the robust linear 
program 

A^A + ^V + CoVW+Cflp < 



SoiX, /i) + M^) + J'ooMS) + ^lol 
COA + DV-7l, + -^o'i</'iW+-^iilp < 










[AsiS)]ijXj + [B5iS)]r,ifij 

is feasible for all 5 G <5 with /i = XliLi /^j' ^0 = 



> 0, 2, j = 1, . . . ,n, i 7^ j 

> 0, z = 1, . . . ,g, j = 1, . . 



'-'0 '-'O 



and 



41 (A, ^) 



A^X + B^Kfi 
ClX + D^Kfi 



A^X + B^Kfi 
A + D'^K^ 



Moreover, in such a case, the controller K is given by 



K 



Ai Vi ••• W?; 



(35) 



(36) 



(37) 



(38) 



and the closed-loop system satisfies H^^iHl^ < 7I |'"^i| U^o- 

Proof: The proof is similar to the one of the results in Section |IV] ■ 
The constraints (|35l) are polynomial in the uncertain parameters S. This is however not the case of the 
constraints (|36] ) that are rational when the system is rational. They can however be easily turned into 
polynomial constraints by finding a common denominator for the left-hand side. Since the sign of the 
common denominator is fixed (otherwise the system would be ill-posed), the rational constraints reduce 
to polynomial constraints on the numerators. 



VIL Solving robust linear programs 

In this section, we address the problem of solving the robust linear optimization problems arising in 
Theorems |2l |4] and |5] A solving scheme based on Handelman's Theorem [|33l is proposed. 



A. Handelman's Theorem 

For completeness, let us first recall Handelman's Theorem: 

Theorem 6 (Handelman's Theorem): Assume iS is a compact polytope in the Euclidean A^-space de- 
fined as 5 := {x G : gi{x) > 0, 2 = 1,...} where the g'i(x)'s are linear forms. Assume also that P 
is a polynomial in variables which is positive on S, then P can be expressed as a linear combination 
with nonnegative coefficients (not all zero) of products of members of {gi}. 

The above theorem states a necessity result regarding the positivity of a polynomial over a compact 
polytope: if it is positive, then we can write it as a linear combination of products of the linear functions 
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defining S. In [|44ll . ||45l it is shown that Handelman's Theorem implies Polya's Theorem [|46l . Sufficiency 
is immediate following similar arguments as for the S-procedure [|22l or sum-of-squares techniques [|47l . 
A linear combination of positive polynomials on S (the products of gi's are positive on S) is indeed a 
positive polynomial on S. Hence, with a suitable choice of the product terms, it is possible to determine 
whether a polynomial is positive over a compact polytope of the Euclidian space. To illustrate the above 
statement, let us consider the following example: 

Example 1: Suppose we would like to characterize all the univariate polynomials p{x) of at most 
degree 2 that are nonnegative on the interval [—1,1]. The basis functions of this interval are given by 
gi{x) = x + 1 and g2ix) = 1 — x. According to Handelman's Theorem, we know that all such polynomials 
write as a linear combination of all possible products gi{xy g2{xy with 1 < i + j <2. Hence, we have 

p[x) = rigi{x) + T2g2ix) + T3gi{x)g2ix) + ngiixf + T5g2{xy ^^^^ 
= X2X^ + Xix + Xo 
where > 0, i = 1, . . . , 5 and 

X2 = n + n- T3, 

Xi = Ti - r2 + 2r4 - 2x5, (40) 

Xo = n+r2 + T3 + T4 + T5. 

Hence, we can conclude on the very general statement that any univariate polynomial of at most degree 2 
that is nonnegative on a bounded interval of the form [a, (3] can be viewed as a point r = (ri, . . . , ts) G M.\. 

B. Equivalent relaxations for robust linear programs 

For simplicity, we will focus here on linear optimization problems depending on a single uncertain 
parameter. The results straightforwardly generalize to the case of multiple uncertain parameters at the 
expense of complex notations. Let us consider the following semi-infinite feasibility problem 

Problem 1: There exists x G M'' such that the inequality 

d 

P(a;,6') := < (41) 

j=0 

holds for all 9 E V := V = {9 E IS. : gi(6) > 0, i = 1, . . . , g[s linear} and where the vectors Pj{x) E 
M^^, j = 0, . . . , c?, are linear in x. 
We then have the following result: 

Theorem 7: The following statements are equivalent: 

i) Problem [T] is feasible. 

ii) There exist an integer 6 > and vectors^ Qk{y) & M^^, k = 1, . . . , 6, such that the finite-dimensional 
linear program 

P,ix) = Zj, j = 0,...,d (Ala) 
Qkiy) < 0, k = 0,...,b (42b) 

^The notation y is here to emphasize that the vectors Qfc's consist exclusively of independent additional decision variables. 
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is feasible in (x, y)eW^x rJ''(^+i) where P{x, 6) = ELo Qk{v)gi{Of, Zj := ^Lo ^^■ik<5fc(?/) and 
the t'jfe's depend on the coefficients of the basis functions giS. 
iii) There exist an integer 6 > and vectors Rk{z) e IR^^, k — 0, . . . , 6 — <i, such that the finite- 
dimensional linear program 





/ 


>o(x)" 












-Ti 






V 






Rb-d{z) 



Rk{z) <0, k 
\ 

< 0. 



0,...,d 



(43a) 



(43b) 



/ 



X 



j,Np(b-d+l) 



where T 



Ti T2 



is feasible in (x, z) e 

T2 e R(d.+i)x{d+i)_ 

Proof: The equivalence between statements i) and ii) is an immediate consequence of Handelman's 
Theorem. The equivalence between statements ii) and iii) follows from simple algebraic manipulations 
allowing to reduce the number of additional decision variables. First remark that the equality constraints 
can be compactly written as 



'Po{x) 




'Qo{y) 




= T 




Pd{x)_ 




_Qb{y)_ 



Note that the matrix T is full-row rank, otherwise it would not be possible to characterize independent 



polynomial coefficients. Using the decomposition T 



Ti T2 



where T2 e 



l)x(d+l) 



nonsingular matrix, the equality constraints can then be solved to get 



'Qi{y) 




/ 


'Po{x) 




Qd+i{y) 










-Ti 




_Qdiy)_ 




V 


Pd{x)_ 




_ Qb{y) _ 



\ 



is w.l.o.g. a 



(44) 



Since the Qi(yys are nonnegative vectors, this is then equivalent to say that the right-hand side of the 
above equality is nonnegative. Finally, posing = Qd+k, k = 1, . . .b — d, we obtain the feasibility 
problem of statement iii). The opposite implication is obtained by reverting the reasoning. The proof is 
complete. ■ 

In the light of the above result, it turns out that it is always possible to represent a robust linear program 
with polynomial dependence on uncertain parameters within a compact polytope as a more complex finite 
dimensional linear program. Handelman's Theorem yields a linear program involving Np{b + 1) + rj 
variables, Np{d+1) equality constraints and Np{b+1) inequality constraints. After reduction, the problem 
has Np{b — d) +77 decision variables and Np{b + 2) inequality constraints. The complexity of the problem 
has then be reduced. The same reasoning applies to the case of multiple uncertainties. 

The question of choosing which products of basis functions to consider is a difficult problem. A brute 
force approach would consider all the possible products up to a certain degree. This is an easy task 
in the case of univariate polynomials. This is however more problematic when considering multivariate 
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polynomials since the number of basis functions b can be very large. A bound on the necessary order for 
the products has been provided in [|44l and generalizes straightforwardly to vector polynomials. 

VIII. Examples 
A. ILCs for several classes of uncertainties 

ILCs for some common operators are discussed below. 

1) Uncertain SISO tranfer function: 

Consider now an uncertain asymptotically stable positive SISO proper transfer function S(s, p) depending 
on constant uncertain parameters p E [0,1]^. By virtue of Corollary [1] the ILC can be expressed as 
0f + (p^Z > where Z G {s(0,p) : p G [0, 

2) Multiplication operator: 

Let us consider in this example the multiplication operator S which multiplies the input signal by a bounded 
and time-varying parameter 6(t) varying arbitrarily within its range of values, i.e. z{t) = T.(w)(t) = 
6{t)w{t). Since the parameter is time-varying, the time-domain version of the ILC must be considered 

r+oo 

/ [^I + flS{9)]w{e)de > 0. (45) 
Jo 

A necessary and sufficient condition is given by ipj + (pJS^t) > for alH > where the scalings cpi and 
ip2 are chosen according to the range of values of 6{t). Note that it is also possible to select polynomial 
scalings verifying (pi{6) = —d(p2{S) in order to saturate the ILC and obtain a better characterization of 
the uncertainty set. 

3) Uncertain infinite dimensional system: 

The same reasoning also applies to general asymptotically- stable positive LTI infinite-dimensional systems 
governed by partial differential equations. Let us consider the heat equation given by 

du 2'^^'" 
dt dx^ 
u{Q,t) = w{t) 

z{t) = u{l,t) 

where u{x, t) is the state of the system, x G [0, 1] the space variable, w{t) the input, z(t) the output and 
a; > an uncertain parameter of the system. The transfer function of the system is given by G{s) = 
Q,g-v^/i^_l_^gv^/tj where a and /3 are real constants determined according to initial and boundary conditions. 
The static-gain G{0) = a + /3 is independent of u and, according to Corollary [1] it is enough to choose 
Lfi = —ip2{a + 13). Stability of interconnections of LTI positive finite-dimensional systems and the heat- 
equation can hence be proved regardless of w G [0, +00). 
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nb. of systems 


{n,p,q) 


I/i-gain 


Loo -gain 


20 


(300,100,150) 


fj. = 12.282, a = 


1.1406 


= 14.186, a = 1.4151 


100 


(50,20,30) 


= 0.53973, cr = 


0.27486 


= 0.50735, a = 0.080446 



TABLE I 

Mean computation time /.t [sec] and standard deviation a [sec] for gain computation 



4) Constant delay operator: 

Let us consider the constant delay operator E(s) = e~^^, h > 0, which is indeed a positive operator. 
The static-gain S(0) is equal to 1 regardless of the value of h. The ILC can hence be easily saturated 
by choosing ipj = —(P2- It is shown in Section IVIII-DI that the consideration of such scalings leads to 
a necessary and sufficient condition for the stability of positive time-delay systems with constant delays, 
recovering then the results of (4E\. 

5) Time-varying delay operator: 

Let us consider now the time-varying delay operator S defined as z{t) = = w{t — h{t)) where 

h{t) > and h{t) < < 1 for all t > 0. In such a case, the operator can only be characterized through 
its Li-gain given by (1 — fj,)^^ under the standard assumption of initial conditions [49|. A suitable ILC 
for this operator is then given by 

POO 

/ tp'^ {w{t) - {1- fi)z{t))dt>0 
Jo 

with (p > 0. Hence ipi = cp and ip2 = —(1 — lAv- 



B. Computation of Norms 

For this example, many linear positive systems have been randomly generated and their induced-norms 
computed on a laptop equipped with an Intel U7300 processor of 1.3GHz with 4GB of RAM. The mean 
computation time and the standard deviation for different systems are gathered in Table IH Note that the 
number of variables is n + 1 and the number of constraints is 2n + p + 1 and 2n + 1 + q for the Li-gain 
and the Loo-gain respectively. Since the number of constraints is larger for the Loo-gain, it is expected 
that it takes longer to compute. It seems important to note that the the induced-norms are different since 
the system is not SISO, so the theoretical analysis in [fT3l does not hold here. 



C. Example 2: Drug distribution 

A frequently used model to analyze the distribution or flow of a drug or a tracer through the human 
body after injection into the bloodstream is given by the following compartmental model [[T]: 



x{t) 



«21 



ai2 
-ai2 



x{t) 



u{t) 



(46) 
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Output matrix C 


T ,A -pain 


T , -pain 


^ r 1 


1 

ail 


1 
ail 


c=[o l] 


■^21 




ail ai2 


aiiai2 


C = dva.g{k\,k2) 


^ll 1^2 a21 
ail a.iiai2 


ail ' aiiai2 J 



TABLE II 

Gains of system (l46l ) for different output matrices 



where Xi, a;2 and u are the compartment corresponding to the blood plasma, the extravascular space (e.g. 
tissue) and the drug/tracer injection. Here, the drug or the tracer is injected directly into the bloodstream 
and the drug is evacuated by the kidneys at rate an > 0. The scalars 012,021 > are the transmission 
coefficients for the drug/tracer between the bloodstream and the extravascular space. Due to the sign 
pattem, the system is positive and asymptotically stable. When the system is SISO the Li-gain coincides 
with the Loo -gain and they are both equal to the static-gain of the system. Different gains for different 
output matrices are given in Table |lll 

D. Theoretical Robustness analysis - Time-delay systems 

We illustrate here the nonconservativeness of Theorem |3] made possible by the saturation of the ILC 
condition. The case of time-delay systems is addressed. 

1 ) Constant time-delay: 

Let us consider the linear positive system with constant time-delay: 

x{t) = Ax{t) + Ahx{t - h) (47) 

for some /i > 0. It is well known that such a system is positive if and only if the matrix A is Metzler 
and the matrix Ah is nonnegative. Rewriting it in an LFT form we get 

x{t) = Ax{t) + AhWo{t) 

zo{t) = x{t) (48) 
wo{t) = Vc{zo){t) 

where Vc is the constant delay operator with Laplace transform Vc(s) = e"'*''. Since the static-gain of 
the delay operator is equal to 1, we can apply Theorem |3] which yields the stability conditions 

X^A + ^J < 0, 

X^Ah + ip^ < 0, (49) 

^J + ip^ = 

that must be feasible for some A G M"^, and (pi,(p2 G K"- They are equivalent to the conditions 

X^A + (^f < and X'^Ah - <Pi < (50) 



21 

which are in turn equivalen1§ to the inequality 

X^{A + Ah)<0 (51) 

well-known to be a necessary and sufficient condition for asymptotic stability of positive time-delay 
systems B^ . [|50l . Note also that this condition is dual to the one derived in [|27l . By solving for A 
instead of (fi, we get the condition 

- ^^A~^Ah < ^1, v^i > (52) 

which is a linear counterpart of the spectral radius condition for the stability of linear time-delay systems, 
see e.g. [|49l . 

2) Time-varying delay: 

Consider now the time-varying delay operator which has been considered in the case of positive systems 
in [|27l . The Li-gain of time-varying operator is (1 — /i)^^, /x < 1 as discussed in Section rVIII-A5l . In 
such a case, the stability conditions are given by 

< »• (53) 

for some A, G Identical conditions can be obtained using a linear Lyapunov-Krasovskii functional 
[|50l of the form V{xt) = X^xi^t) + V5^//_;j(^-) x{s)ds with A, (/? > 0. Note that in the robust formulation, 
(1531 ) is feasible only if (/? > 0. Note also that the conditions can be merged into the novel single inequality 

X^iil- fi) A + Ah) <Q (54) 

where we can see that the time-varying delay penalizes the stability by scaling down the matrix A. The 
time-invariant case is retrieved for /i = 0. 

While the above results are based on the Li-gain, it seems interesting to analyze the stability using 
the Loo-gain. The Loo-gain of the time-varying delay operator is equal to 1 [tSTI . We may then use a 
small-gain argument to prove asymptotic stability of the time-delay system by applying a scaled version 
of Lemma m Hence, the time-delay system is asymptotically stable provided that there exist A, G M"_|_ 
such that the conditions 

^A + Ahip < and A - v?^ < (55) 

hold. They are equivalent to the condition (^4 + Ah)X < (stability for zero delay) or the condition 
—A^^Ahip < ip (small-gain condition) which is the Loo-gain counterpart of (|52] ). It is important to precise 
that the stability condition is identical to the one for systems with constant delay. This result is then 
much stronger than the Li -based one since the delay-derivative bound does not have any negative impact 
on the stability. As discussed in [ISTTl . while results using the Li-norm relate to Lyapunov-Krasovskii 
functional, those based on the Loo -norm connect to Lyapunov-Razumikhin functions. Note however that 
no Lyapunov-Razumikhin results for positive systems have been reported so far. 

^To see the equivalence, note that iSOl imphes i51t by summation. To prove the converse, assume ( I51t holds and defining (pj — 
Ah — {A + Ah) 12 > in J50t makes the conditions negative. 
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E. Numerical robustness analysis - A systems biology example 

Let us consider the example of a simple gene expression process described by the following model 



u{t) 



(56) 



Xrit) — 7r Xrit) ^ 1 

Xp{t)\ [ kp -7pJ \xp{t)\ [O 

where > is the mean number of mRNA in the cell, > is the mean number of protein of interest 
in the cell. Above u{t) >Q is the transcription rate of DNA into mRNA, 7^ > is the degradation rate 
of mRNA, fcp > is the translation rate of mRNA into protein and 7^ > is the degradation rate of the 
protein. The parameters are assumed to be uncertain and given by 7r = 7° + Si 7^, kp = + £2 k^ and 
7p = 7p + £3 7p, where £j G [— 1, 1], z = 1, 2, 3. We are interested in analyzing the Loo-gain of the transfer 
from u Xp over the set of all possible systems. In other terms, we would like to analyze the impact of 
the maximal value of the transcription rate to the maximal value of the mean number of proteins. 
The overall system can be rewritten as 



Xrit) 


= Auie) 


Xrit) 


+ 


1 


Xpit) 




Xpit) 








uit) 



(57) 



where e = col(£i, £2, ^3) and 
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0" 




"0 


0" 
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-7° 






+ £2 


k^ 




+ £3 




> 















Auie) 



Note that the theoretical computation of the Loo-gain is rather difficult since the static-gain depends 
rationally and in a nonconvex way on the uncertain parameter vector e. However, by exploiting the fact 
that the set := {Auie) : e E [—1,1]^} is convex, the conditions or equivalently (fTOl ). can be checked 
without using the LFT formulatiorQ although in this case both methods are equivalent provided that the 
scalings are chosen accordingly. Checking indeed the conditions ^ over the whole set £/ is equivalent to 
checking them over the set of its vertices, i.e. ^/^ := {A^ie) : £ G { — 1, 1}^}- This follows from standard 
convexity argument. 

Let us consider 7° = 1, Ai^ = 2 and 7° = 1 for numerical application. We also assume that the 
parameters are known up to a percentage N e [0, 1) of their nominal value, hence 7^ = A^7°, 7^ = A^7p 
and kp = Nk^. The results are summarized in Table In] where we can see that the proposed approach is 
rather accurate. 



F. Numerical robustness analysis - A polynomial system example 

Let us consider the uncertain system with constant parametric uncertainty 6 E [0, 1]: 

x{t) = iA'> + 6A' + 6^A^)xit) + iE'> + 6E' + 6^E^)wiit) 
z^it) = iC^ + + S^C^)xit) + iF'> + + 6^F^)wiit) 



(58) 



Remember that the interest of the LFT is to convexify the problem. 
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N 


Loo -gain 


Theoretical 





2 


2 


0.1 


2.7162 


2.7161 


0.3 


5.3063 


5.3062 


0.5 


12.0003 


12.0000 


0.7 


37.7783 


37.7779 



TABLE III 

Computed and theoretical Loo-gains of the transfer u ^ XpOF the uncertain system dSTl i 



MS) 


<fi2{5) 


constraints 


computed Li-gain 


time [sec] 




V2 + H^\5 


v\ = -'^2, fi = -vh 


133.95 
133.95 
94.167 


2.7844 
3.829 
4.2758 



TABLE IV 

Li-GAIN COMPUTATION OF THE TRANSFER Wi Zi USING THEOREm[2]- EXACT Li-GAIN: 92. 835 ^ 



with the matrices 
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The system can be rewritten in the LFT form (fTTl with matrices = A^, Eq 



E, = E\ Ci = C\ F, 



10 — 

In ' 

On 



(Jl (J2 pi p2 



Fu 



F 



00 



F° and 



Onxp Onxp 
Onxp Ojjxp 



Opxn Opxn 
Opxn Opxn 



0. 







^01 



The matrices of the LFT form (|32] ) are given by Foo = Fqo, Coo 



pxp 

Co, Fqi 



Foi, Eq 



""nxp 
"I 

Op _ 



(59) 



A^ E^ E^ 



A 



2T 



(60) 



C^^ c 



and Fi 



10 



E^' E 



2T 



:^2T 



Using Theorem |2] and Theorem 



4] (together with the complexity reduction technique of Section IVIII) 
and different forms for the scalings ipi and ip2, we obtain the results gathered in Tables |IV] and |Vl We 
can see that the Li-gain is not very well estimated compared to the Loo -gain for parameter independent 



24 







constraints 


computed _/ji-g3.in 


time [sec] 


^? 

As 


Ai 
ip2 + 


V'? > 0, v'? + (^2 > 

V^l = -9=2 


86.195 
86.195 
82.025 


0.68989 
1.4629 
1.7509 



TABLE V 

Loo-GAIN COMPUTATION OF THE TRANSFER Wi -> Zi USING THEOREm|4]- EXACT Lcxj-GAIN: 82.0249 

scalings. Using scalings of degree two considerably reduces the conservatism of the approach. Indeed, 
we are able to estimate accurately the Loo-gain while little conservatism persists for the Li-gain. It is 
important to note that the above numerical results also hold in the case of time-varying parameters (see 
Section IVIII-A2I) . 

IX. Conclusion 

By relying on the fact that the Li-gain and Loo-gain are related to the static-gain matrix of the positive 
system, a linear programming approach has been proposed to compute them under the assumption that the 
input signals are nonnegative. Interestingly, the results turn out to be also valid for inputs and states with no 
definite sign. Based on this result, stabilization conditions with performance constraints have been obtained 
for different classes of controllers. Using then the newly introduced Integral Linear Constraints, robust 
stability analysis results are provided together with generic robustness results, again relying on static-gains 
of positive systems. It is shown that the results are nonconservative in the case of UTI positive uncertainties 
with fixed static-gain matrix. Stabilization results are finally obtained and expressed as (robust) linear 
programming problems. An exact relaxation scheme based on Handelman's Theorem is proposed in order 
to obtain finite-dimensional feasibility problems. Several examples illustrate the approach. 
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